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Abstract. The wave equation utt — c^Uxx is generally regarded as a linear approximation 
to the equation describing the amplitude of a transversely vibrating elastic string in the 
plane. But, as is shown in 0, the assumption of transverse vibration in fact implies that 
the wave equation describes the vibration precisely, with no need for approximation. We 
give a simplified proof of this result, and we generalize to the case of an elastic string 
vibrating (transversely or not) in a Riemannian surface M. In the more general setting, 
the assumption of transverse vibration is replaced by the assumption of "perfect elasticity," 
and we show that the wave map equation Vu^Ut = c^Vu^Uj, gives a precise description of 
the vibration of a perfectly elastic string in M, with no need for approximation. Finally, we 
give examples describing the motion of various vibrating strings in E^, 5^, and H^. 



1. Introduction 

One of the most familiar partial differential equations in all of mathematics is the wave 
equation for a vibrating string, 

(1-1) Utt = C^U^x, 

which describes the transverse vibration of an elastic string the in the plane. The equilibrium 
position of the string is assumed to be an interval [0, L] along the x-axis, and the function 
t) represents the vertical displacement at time t of the point on the string corresponding 
to the point with equilibrium position (x,0). 



A wide variety of methods have been used to derive equation ( 1.1 ) from physical principles; 
we will critically review some of the best-known derivations ([7], [6], [12]) in ^ The methods 
discussed here vary in the specific techniques used during the derivation process, but they all 
make use of simplifying assumptions and/or approximations (e.g., constant tension, small 



amplitude vibrations), with the result that equation (1.1) is deemed to be only a linear 
approximation to the "true" wave equation. Exceptions to this approach in the literature 
are rare; a few can be found in P, [2], [8], [T6] . 

In ||3} we will start from scratch and carry out a derivation of the wave equation from 
physical principles, without any of the simplifying assumptions made in most derivations. Our 
derivation is along the lines of that given in [2] , but it is more straightforward. Remarkably, 



it turns out that all the nonlinearities cancel each other out, and equation (1.1) is the 
correct equation for the motion of a transversely vibrating elastic string — not "just" a linear 
approximation. In light of the simplicity of our derivation, we find it remarkable that this 
fact is not better known than it appears to be, and that the presumed "nonlinear string 
equation" has persisted for many decades in textbooks. 



In ^ we will consider generalizations of equation (1.1 ) to the case of an elastic string with 



not-necessarily-transverse vibration, both in the Euclidean plane and in an arbitrary 
Riemannian surface M, where the notion of "transverse" vibration is not well-defined in 
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general. In this more general setting, the assumption of transverse vibration is replaced 



by the assumption of "perfect elasticity" (cf. Definition 3.1). We will see that with this 
assumption, the motion is governed by the wave map equation 

(1.2) Vu.ut = c'Vu^u, 

for the function u : [0, L] x M — )• M, where V represents the Levi-Civita connection of the 
Riemannian metric on M. (We refer the reader to [13], [IT], and [H] for information on wave 
maps.) Again, while the wave map is often portrayed as a "first approximation" (see, e.g., 
[13] ) to the motion of an elastic string, we will show that in fact no simplifying assumptions 



are needed in order to derive the wave map equation (1.2) from the physical principles 
governing the motion of the string. In ^ we give some examples of solutions to (1.2) in 
non-flat Riemannian surfaces. Note that here the partial differential equation governing the 
motion of the string on a non-fiat Riemannian surface is intrinsically nonlinear. This arises 
from the non-Euclidean nature of the ambient space in which the vibrations are taking place. 

2. Many derivations, many assumptions 

Consider an elastic string of constant linear density p, whose equilibrium position lies along 
the interval [0, L] in the x-axis. Suppose that the string is allowed to vibrate only in the 
transverse direction, so that at time t, the position in the xy-plane of the point on the string 
corresponding to the point with equilibrium position (x, 0) is given by {x,y) = {x,u{x,t)). 

2.1. The Courant-Hilbert derivation. In Courant and Hilbert's classic text [7], the wave 
equation is derived from Lagrangian mechanics and the calculus of variations. Simplifying 
assumptions include: 

• The magnitude T of the tension in the string is assumed to be constant. 

• The vibrations are assumed to have "small" amplitude, so that the quantity >/l + 
can be approximated by its first-order Taylor polynomial 1 + ^m^. 

The kinetic energy K of the string is given by 

K = p ul dx, 



JO 

while the potential energy P is assumed to be proportional to the increase in the string's 
length compared to its length at rest; i.e., 

P = t(^J^ ^^lT^^dx-l}j uldx. 

This leads to the Lagrangian functional 

C = K -P= ! [pul - Tul) dx, 
Jo 

whose Euler-Lagrange equation is 

putt - Tu^^ = 0. 



Setting c = - yields the wave equation (1.1) 
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2.2. The Coulson-Jeffrey derivation. Coulson and Jeffrey's well-known text [B] takes 
a different approach that avoids the calculus of variations, and instead considers the forces 
acting on a small segment of the string, corresponding to the interval [x, x + Ax]. Simplifying 
assumptions include: 

• The magnitude T of the tension in the string is assumed to be constant. 

• While it is noted that the condition for T to be approximately constant is that "the 
wave disturbance is not too large," this assumption is not made explicit as in the 
Courant-Hilbert derivation until fairly late in the process. 

The equation of motion for this segment of the string is given by Newton's law of motion 
F = ma; as the vibration is assumed to be transverse, only the vertical component of the 
motion is considered. 

For each fixed t, let d{x,t) denote the angle between the tangent vector to the string at 



the point {x,u{x,t)) and the horizontal, as shown in Figure 2.1 



T{x + Ax,t) 




X X + Ax 

Figure 2.1. String vibrating transversely in the plane 

The vertical component of the force exerted on the right-hand end of the string segment is 
given by Tsin6'(x-|- Ax, t), while the vertical component of the force exerted on the left-hand 
end of the string segment is given by — Tsin6'(x,t), so that the total vertical force acting on 
the string segment is 

F = Tsine(x + Ax,t) -rsin^(x,t). 
The arc length of the string segment is approximately 

As = a/1 + ul Ax, 

and so the mass of the string segment is assumed to be 

m = pAs = pa/1 + ul Ax. 

Note that this assumption is incorrect: the mass of this string segment at rest is equal to pAx, 
not pAs, and conservation of mass implies that the mass remains equal to pAx regardless of 
how the string is stretched. (This error appears frequently — although not universally — in the 
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literature; it is generally inadvertently remedied at a later stage in the derivation where u'^ 
is assumed to be small, and therefore As ^ Ax.) Since a = Uu, Newton's equation becomes: 

putt^/l + ul Ax = Tsm9{x + Ax,t) — T sin 9{x,t). 
Dividing by Ax and taking the limit as Ax — ?■ yields 

„ d 

putt\/l + ui 



(2.1) 



T-(sin^(x,t)) 
T(cos^(x,t))0x- 



The trigonometric identities 



(2.2) 

imply that 

Therefore, 
(2.3) 



tan 6' 



sin^ 



d_ 
dx 



(tan^) 



sec 



cos 6 



ut 



ut 



and equation (2.1) is equivalent to the non/mear equation 

TUxx 



(2.4) 



PUtt 



(If the mass of the string segment had been correctly represented as pAx instead of pAs, 
equation (2.4) would instead have become 

2.5 pUtt = . 2^3/2 - 

Only at this stage in the derivation do Coulson and Jeffrey apply the assumption that is 
small, so that equation (2.4) may be approximated by the linearized form (1.1), with = ^. 

2.3. The undergraduate version. Yet another variant on this derivation appears in many 
undergraduate PDE textbooks; as an example, we present a slight variation on the derivation 
given in [12] . In this version, assumptions about the tension T and the magnitude of Ux are 
postponed until later in the process. As in the Coulson- Jeffrey derivation, this argument 
is based on Newton's equation applied to the string segment corresponding to the interval 
[x,x + Ax], but this time Newton's equation appears in its vector form F = ma. With the 
same notation as in ^2.2 (and using m = pAx), the vertical component of Newton's equation 
becomes 

(2.6) puttAx = T(x + Ax,t) sin^^(x + Ax,t) — T{x,t) sm9{x,t), 
while the horizontal component becomes (due to the assumption of transversality) 

(2.7) T(x + Ax, t) cos e{x + Ax, t) - T{x, i) cos ^(x, t) = 0. 

First consider equation (2.7). Dividing by Ax and taking the limit as Ax — ?■ yields 
d 



dx 



(T(x,t) cos^(x,t)) = T^(cos^) -T(sin^)^^. = 0. 
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Now the simplifying assumption is made that 6 is "small," and therefore 



sin6' ^ 0, 



cos 9 



hence, ~ 



0. So T is now assumed to be approximately constant, but in this case the 
assumption of constant tension is a consequence of the hypothesis of transverse motion and 
the assumption that the vibrations are "small," rather than a separate assumption. 



With these assumptions in hand, equation (2.6) becomes 

putt^x = Tsin6{x + Ax,t) — T sm6{x,t), 



and essentially the same argument as in 32. 2 1 shows that, with the assumption that ~ 0, 



this equation is approximately equivalent to equation (1.1). 



3. A MORE ACCURATE DERIVATION 



In this section we present an alternative derivation for equation (1.1). We will follow 



the general strategy of the derivation of ^2^, but we will assume only that the vibration is 
transverse, so that motion occurs only in the vertical direction. No a priori assumptions will 
be made about either the function T{x, t) (except that it is assumed to be differentiable) or 
the magnitude of the vibrations. (Derivations using only these assumptions are given in [2] 
and [8]; ours is similar, but more straightforward.) 

As in ^2.3, we apply Newton's equation F = ma to the string segment corresponding to 



the interval [x,x + Ax], which leads to equations (2.6) and (2.7). Dividing both of these 
equations by Ax and taking the limit as Ax — ?■ yields the equations 

d 

(3.1) pwit = — (T(x,t) sm6'(x,t)) , 



(3.2) 



d_ 
dx 



dx 

(T(x, t) cos 6'(x, t)) 



First consider equation (3.2). Write this equation as 

T^{cos9) -T{sme)9^ 



0. 



0, 



and substitute in the trigonometric identities ( |2.2 ), (2.3) to obtain 

1 , '^-^ x'^^xx 



(3.3) 



T 



ut 



(1 + «2)3/2 



0. 



Equation (3.3) may be regarded as a separable differential differential equation for the func- 



tion T{x,t): rewrite equation (3.3) as 



T 

Y 



UxUx 



Now integrate and exponentiate to obtain: 

(3.4) T(x,t) = C{t)^l + ul, 

where C{t) is an arbitrary function of t. 



Observe that, for each fixed t, equation (3.4) says that the magnitude of the tension of 



the string at the point (x, ^(x, t)) is proportional to the derivative of the arc length function 

fX 

s{x,t)= / ^/l + Ux{x,t) dx 
Jo 
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with respect to the curve parameter x. Assuming that the physical parameters of the string 
are not changing in time, it is reasonable to assume that the "constant" of proportionality 
C{t) is, in fact, equal to a constant Tq, which represents the tension of the string in its 
equilibrium position. The quantity 

ds n 



dx 



ut 



may be regarded as the "stretching factor" of the string at the point (x, u{x, t)), and, following 
we make the following definition: 



Definition 3.1. A string is called perfectly elastic provided that, when the string at tension 
To is stretched by a factor of a, the tension becomes equal to T = Tqct. 



Thus, equation (3.4) may be interpreted as saying that the hypothesis of transverse vibra- 



tions implies that the string is perfectly elastic, merely as a consequence of Newton's law. 
This observation is significant enough that we state it as: 

Proposition 3.2. The tension of an elastic string vibrating transversely in the plane must 



satisfy the perfect elasticity condition (3.4). 



Remark 3.3. Perfect elasticity is a special — and idealized — case of the more familiar notion 
of "linear elasticity." Linear elasticity refers to a string that satisfies Hooke's Law, which 
in its infinitesimal version says that, at a point where the string's stretch factor is equal 
to cr, the tension in the string is proportional to {a — ctq), where ctq is the value of ^ for 
which the tension is equal to zero. Perfect elasticity is the special case where ctq = 0, which 
mathematically corresponds to a string which shrinks to a point in the absence of external 
forces. As a practical matter, this condition is almost certainly never precisely satisfied 
by real physical strings; however, a rubber band or a spring stretched to several times its 
natural length may satisfy it to a reasonable approximation. (Physical data supporting this 
assertion may be found in [S].) Consequently, the vibrating motion of real physical springs 
can never be truly transverse, but only approximately so. 



Continuing with the derivation, substitute the trigonometric identity (2.2) and the expres 



putt 



sion (3.4) for T(x,t) into equation ( 3.1[ ) to obtain 

dx 

d_ 
dx 

C{t)u, 



(3.5) 




With the assumption that C{t) = Tq, equation (3.5) is equivalent to (1.1), with = ^. 

4. Vibrating strings in Riemannian surfaces 

Now consider the problem of an elastic string of constant linear density p vibrating in 
a Riemannian surface M. Even the innocuous-sounding assumption of "constant linear 
density" merits some examination. If a string is represented by a parametrized curve 

u : / ^ M, 



6 



for some interval / C M, then the density of the string is represented by a 1-form dfi on /, 
defined by the condition that the mass of any segment u([a, b]) of the string is given by 

m(a, b) = I dfi. 

J a 

"Constant hnear density" is generally taken to mean that d^ = pdx for some constant p, 
but in fact this condition is dependent on the parametrization of the string. A string with 
any smooth, nonvanishing density 1-form dp can always be given a reparametrizion u(m) 
with the property that dp = pdm for some constant p > 0. (In fact, by taking m to be the 
function 

r 

m{x) = / dp, 



JO 

we can arrange that p = 1; however, it is often convenient to allow other values of p, so we 
will allow p to be any positive constant.) 

Definition 4.1. Let I G M. be a connected interval, and let m be a local coordinate on I. A 
parametrized curve u : I ^ M in a Riemannian surface M with a specified density 1-form 
dp on I is said to be parametrized by constant density if dp = p dm for some constant p > 0. 

As observed above, a string with any smooth, nonvanishing density 1-form can be para- 
metrized by density, so the assumption of "constant linear density" must mean something 
more than simply having the string parametrized by constant density. In the classical case 
of a string vibrating transversely in the plane, the condition of constant linear density is 
characterized by the fact that the equilibrium position of the string has a parametrization 
which is simultaneously a parametrization by constant density and a parametrization by 
arc length. As the vibrating string expands and contracts, however, the parametrization 
{x,u{x,t)) for the position of the string at a given time t generally does not remain an arc 
length parametrization for the curve — but it does remain a constant density parametrization. 

Next, consider the notion of "transverse vibrations." In the classical case of the string 
vibrating in the plane, "transverse" is taken to mean "in the vertical direction." This notion 
relies crucially on the canonical local coordinates (x, y) on the plane, and on the vertical 
vector field defined on any open neighborhood in the plane. But for a general Riemannian 
surface, there are no canonical local coordinates, and hence there is no obvious way to define 
the notion of a "transverse" vibration. (Even if we attempted to accomplish this by defining 
a "transverse" vector field orthogonal to a curve representing the equilibrium position of 
some particular string, there is still no canonical way to extend such a vector field to an 
open neighborhood of the given curve segment.) 



Therefore, in light of Proposition 3.2 we replace the assumption of "transverse vibration" 



with an assumption of perfect elasticity, which does make sense even in the more general 



setting of a Riemannian surface. But this too requires some care: according to Definition 3.1 



a string is perfectly elastic if the tension in the string is directly proportional to the stretch 
factor jf-. But like the constant linear density condition, this condition is dependent on the 



parametrization of the string. Since Definition 3.1 was formulated under the assumption of 
a constant density parametrization for the string, we now make the following refinement in 
order to make this assumption explicit: 



Definition 4.2. A string in a Riemannian surface M with a constant density parametriza 

ds 
dm 



tion u(m, t) is called perfectly elastic for a = in the range a < a <b provided that, when 
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the string at tension Tq is stretched by a factor of a E [a,b], the tension becomes equal to 
T = Toa. 

In order to derive the equation of motion for the general case of a string vibrating in a 
Riemannian surface M, we will assume that: 

• the string is parametrized by constant density; 

• the string is perfectly elastic; 

• the motion of the string is governed by Newton's Law. 
Corollary |4.4 will show that these conditions suffice to guarantee that a constant density 



parametrization for a string in equilibrium position is also a constant speed parametrization 
with respect to an arc length parameter along the curve segment. This is the natural analog 
to the "constant linear density" assumption in the classical case; the fact that it follows from 
the assumption of perfect elasticity suggests that perhaps the notion of perfect elasticity is 
more subtle than it might appear. 

We are now ready to state our main result: 

Theorem 4.3. Let M be a Riemannian surface with Levi-Civita connection V, and let 
u : [0, L] X M — )■ M 5e a smooth map, where for each fixed t E R, the curve u([0, L] x {t}) 
represents the position of a string vibrating in M at time t. Suppose that the string is 
parametrized by constant density and perfectly elastic, and that its motion is governed by 
Newton's Law. Then the map u(m, t) satisfies the wave map equation 

(4.1) Vu.Ut = C^Vurr^Um- 

Corollary 4.4. The equilibrium position of a string in a Riemannian surface M satisfying 
the conditions of Theorem 4-3\ consists of a segment of a geodesic curve in M . Furthermore, 



the given constant density parametrization of this curve segment is also a constant-speed 
parametrization with respect to an arc length parameter along the curve segment. 



Proof of Corollary 4-4' Any equilibrium solution is independent of t; therefore VutU^ = 0. 



Then (4.1) implies that the solution curve u(m,t) = u(m) satisfies the geodesic equation 



Vu„u^ = 0, 



which implies that u(m) is a constant-speed geodesic segment in M. □ 



For the sake of clarity, we will prove Theorem |4.3| in two steps: first we give a proof in the 
case that M is the fiat plane E^, where we can take advantage of the canonical identification 
between tangent planes at each point of E^. Then we will indicate how the proof may be 
generalized for an arbitrary Riemannian surface M. 



Proof of Theorem 4-3. First suppose that M = E^, and let 

u(m, t) = (x(m, t),y{m, t)) . 

For fixed t, consider the string segment corresponding to the interval [m, m + Am] and the 
(vector-valued) tension forces T(m,t),T(m -|- Am,t) acting at the endpoints of the string 



segment. (See Figure 4.1) According to the assumption of perfect elasticity, the magnitude 



of the force vector T(m, t) is equal to 

T{m,t) = To —{m,t) = Tq (^a/x^ + 



{m,t) 
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Figure 4.1. Perfectly elastic string vibrating in the plane 



for some positive constant Tq; similarly, the magnitude of the force vector T(m + Am, t) is 
equal to 

T{m + Am, t) = Tq ^{m + Am, t) = Tq [V^l + vl] 
The direction of T(m + Am, t) is given by the unit tangent vector 



(m+Am,t) 



t{m + Am,t) 



(m+Am,t) 

to the curve at the point u(m + Am, t), while the direction of T(m, t) at the point u(m,t) 
is given by 



— t(m, t) 

Therefore, we have 

T(m,i) = -T{m,t)t{m,t) = 
= -TQUm{m,t), 



-1 



- (^0 V^^ + yl) 



{m,t) 
1 



{•^mi Um) 



T(m + Am, t) = T(m + Am, t)t{m + Am, = (Tq y/xl^+yl) \ 



= To Um(m + Am, t). 
Thus the total force acting on the string segment is 

(4.2) F = T(m, t) + T(m + Am, t) = Tq (u^(m + Am, t) - u^(m, t)) . 



(m,t) 



(m+Am,t) 
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Since u(m, t) is assumed to be a constant density parametrization, the mass of the string 
segment is given by p Am. Now Newton's equation becomes 

(4.3) (p Am) utt = To [um{m + Am, t) - u„(m, t)] . 

Dividing by Am and taking the hmit as Am — i- yields 

p Utt = To 

Taking = — yields 



p 

Uft = VLr 



which is equivalent to (4.1) for the standard flat connection V on E . 

For a general Riemannian surface M, we must make the following modifications to this 
proof: 

• In Newton's equation, the acceleration vector utt must be replaced by VutU^. 

• The expression 

T(m,t) + T(m + Am,t) 



in (4.2) no longer makes sense because the tangent vectors T(m,t), T(m + Am,t) 
are based at different points of M. In order to remedy this, we must replace the 
vector T(m + Am, t) by its parallel transport backwards along the curve u(-, t) from 
u(m + Am, t) to u(m, t). 



Am — 7- yields the wave map equation (4.1) 



Substituting these changes into equation (|4.3|), dividing by Am and taking the limit as 

□ 

5. Examples 



In this section we compute several examples of solutions to equation (4.1) (with c = 1) in 
various Riemannian surfaces. These examples illustrate how solutions may be affected not 
only by the geometry of the surface M, but also by the choice of a density 1-form dp, along 
the string. (The choice of dp, is implicit in the choice of parametrization for the initial data 
curve, as this curve is assumed to be parametrized by constant density.) 

For each example we will specify a Riemannian surface M and an initial curve segment 



u(m, 0) = Uo(m) in M. We will give an explicit representation for the PDF system (4.1) 
in terms of local coordinates on M, and we will numerically solve this system for the given 
initial data, assuming fixed endpoints and zero initial velocity]^ 

Example 5.1. Let M = be the flat plane. With local coordinates (x, y) and the standard 
flat metric on E^, the wave map equation (4.1) is equivalent to the system 



(5.1) Xtt — Xmmi Utt — ymm 

for the functions x(m, t),y{m, t). 

Suppose that the initial curve is the graph of y = sin(7rx) for < x < 1. The classical 
case of "constant linear density" along this curve is represented by the parametrization 

x{m, 0) = m, y{m, 0) = sin(7rm), < m < 1. 



"'^AU numerical solutions were computed in Maple 16 using the "pdsolve/numeric" function with standard 
options. 
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This parametrization leads to the classic solution 

x{m, t) = m, yiiTL, t) = cos(7rt) sin(7rm). 



Some curves in this evolution are shown in Figure 5.1 the initial curve is drawn as a thick 
curve, while subsequent curves in the evolution are thinner. 




Figure 5.1. Evolution of "constant linear density" string in 

Now suppose that instead of the density 1-form = dx, we choose dfi = 2{x + ^) dx, so 
that the density of the string is an affine linear function of the initial x-coordinate of the 
string. This corresponds to choosing m = {x + 

x{m, 0) = ^/Tn — ^, yijn, 0) 



^^"^j which yields the parametrization 

\)), \<m<l 



2 

sinfvr 



for the initial data curve. The system (5.1) can be solved explicitly, either by Fourier series 



or by D'Alembert's formula; some curves in this evolution are shown in Figure 5.2 Note 




Figure 5.2. Evolution of string with density d^i = 2[x + \) dx in 

that, regardless of the initial curve, the general theory of the linear wave equation guarantees 
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that the fixed-endpoint problem will have periodic solutions. The period, however, depends 
on the parametrization: if the initial curve has parametrization 

x(m, 0) = Xo(m), y{m,0) = yo{m), a<m<b, 

then the vibration of the string will have period equal to 2(6 — a). 

Example 5.2. Let M = S"^ be the unit sphere in M.^, with local parametrization 

TT TT 

(x, —)■ (cos X COS sin x cosy, sin y), < x < 27r, ^2^^^ 2' 
The standard metric on S'^ is given by 

ds"^ = (cos^ y) dx^ + dy'^, 
and the wave map equation on S"^ is equivalent to the system 

(5.2) Xtt = Xmm ~ 2tan?/ {XmUm ~ ^tUt), 

ytt = Vmrn + sin y cos y (x^ - xl). 

This is a nonlinear system, and its solutions are not periodic in general. 

By analogy with our previous example, we will consider "sinusoidal" initial curves. First 
suppose that the initial curve is described in local coordinates on S*^ by the equation 

1 /27r \ 3 
?/ = - sm —X , < X < 



2 V 3 / ' - - 2' 

with parametrization 

N 1 /^27r \ 3 
x(m, 0) = m, y[fn^ 0) = - sin I — m I , < m < -. 

2 Y 3 y 2 

Figure |5.3 shows the evolution of the curve through three cycles in the vibration; as in the 
last example, the thick curve in each picture represents the initial curve for that portion 
of the evolution. Note that this solution is only quasi-periodic, as the shape of the wave 
changes with each successive cycle. 

The effect of the nonlinearity becomes more significant if we increase the amplitude of the 
initial curve: suppose that the initial curve is described in local coordinates on S'^ by the 
equation 

/27r \ 3 
1/ = sm I yx 1 , < X < -, 

with parametrization 

x(m, 0) = m, yiTTi, 0) = sin i "^^"^ j ' < m < ^. 



Figure |5.4| shows the evolution of the curve through three cycles in the vibration; note that 
the shape of the wave varies dramatically and even becomes non-embedded at times. 

Example 5.3. Let M = be the Poincare upper half plane 

= {(x,y) gm^ I 2/>o}, 

with its standard metric 

, ^ dx'^ + dy'^ 
ds^ = . 

y 
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Figure 5.4. Evolution of large amplitude string in S"^ 
The wave map equation on is equivalent to the system 

(5.3) 



ytt = Umm + 



y 

{xl 1^ xj) - {yl - yl) 



As in the previous example, this is a nonhnear system, and its solutions are not periodic in 
general. 

We will consider several different initial curves. First consider the horizontal line segment 
y = 1 for < X < 2, with parametrization 

x(m, 0) = m, yiiTL, 0) = 1, < m < 2. 

(Recall that this curve is not a geodesic in EP.) Figure 5.5 shows the evolution of the curve 
through three cycles in the vibration; the solution appears almost (but not quite) periodic. 





Figure 5.5. Evolution of horizontal string in EP 

Next, consider a "sinusoidal" perturbation of a vertical geodesic segment. In order that 
the initial curve appear approximately sinusoidal with respect to the Poincare metric, we 
will take the curve ^ 

x=-ysm{y-l), l<Z/<vr + l, 

as our initial curve, so that the horizontal line segment from the point {x{y),y) to the point 
(0,1/) has length ^sm{y — 1). (Note that this is still only approximately sinusoidal, since 
the horizontal line segment joining these two points is not a geodesic.) Furthermore, we 
will parametrize this curve according to an arc length parametrization for the corresponding 
vertical line segment along the y-axis: 

a;(m,0) = Je™sin(e'" - 1), y{m,0)=e"', < m < ln(7r + 1). 



Figure |5.6| shows the evolution of the curve through three cycles in the vibration; one in- 
triguing feature of this evolution is that within each cycle, the string moves to the left much 
more rapidly than to the right. 
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Figure 5.6. Evolution of small amplitude string in 



As in the previous example, the effect of the nonlinearity becomes more significant if we 
increase the amplitude of the initial curve: suppose that the initial curve is described by the 
equation 

X = ysm{y — 1), 1 < y < tt + 1, 

with parametrization 

x(m,0) = e"'sin(e™ - 1), ?/(m,0)=e'", < m < ln(7r + 1). 
Figure |5.7 shows the evolution of the curve through three cycles in the vibration; now the 



shape of the wave varies dramatically and even becomes non-embedded at times. 

6. Conclusion 

In this paper we have shown that the classical, frictionless, elastic string under tension 
undergoing transverse vibrations in the Euclidean plane is exactly modeled by the linear 
wave equation, 

(6.1) Uu = C^Uxx, 

regardless of the size of the vibrations, and that this PDF is not a small amplitude approxi- 
mation of a "nonlinear string equation." More generally, if the arena of the vibration changes 
from the Euclidean plane to an arbitrary Riemannian surface, such the 2-sphere S"^ or the 
Poincare half-plane H^, then the motion of the string is governed by the partial differential 
equation for wave maps into Riemannian surfaces: 

(6.2) Vu.ut = cVu.u,. 

The wave map equation is intrinsically nonlinear whenever the Gauss curvature of the target 



space of the wave map is non-zero. Equation (6.1) is essentially the wave map equation 
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Figure 5.7. Evolution of large amplitude string in 



for wave maps into the Euclidean plane. Equation (6.2) is a source of fascination for both 



the mathematician and physicist. If forces cease to act on the string, then the left hand 



side of (6.2) will vanish, and the resulting differential equation will be that of a geodesic 



of the Riemannian surface. Hence (6.2) generalizes the notion of geodesies and provides a 
very interesting class of hyperbolic partial differential equations for geometric analysis ([11]. 
[H]). For the mathematician, there are still many basic questions that remain open. On the 
other hand, the same partial differential equation arises, for instance, in elementary particle 
physics ([0]) and general relativity {[IQl, ^). 

One natural question to consider is which Riemannian surfaces give rise to wave maps with 
particularly nice properties. In p] we have begun an exploration of wave maps for which 
the corresponding wave map system is Darhoux integrable; we plan to delve more deeply 
into this topic in future work. Another area which we hope to explore in the near future is 
that of sub- Riemannian geometry; the context in which we have discussed wave maps here 
as parametrized curves evolving in time still makes sense in sub-Riemannian geometry, and 
questions such as global existence of solutions and integrability are of interest. 
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